Tarea 3

Preparativos

primates <-
  st_read(
    "https://raw.githubusercontent.com/gf0604-procesamientodatosgeograficos/2021i-datos/main/gbif/primates-cr-registros.csv",
    options = c(
      "X_POSSIBLE_NAMES=decimalLongitude",
      "Y_POSSIBLE_NAMES=decimalLatitude"
    ),
    quiet = TRUE
  )

st_crs(primates) = 4326


cantones <-
  st_read(
    "https://raw.githubusercontent.com/gf0604-procesamientodatosgeograficos/2021i-datos/main/ign/delimitacion-territorial-administrativa/cr_cantones_simp_wgs84.geojson",
    quiet = TRUE
  )
provincias <-
  st_read(
    "https://raw.githubusercontent.com/gf0604-procesamientodatosgeograficos/2021i-datos/main/ign/delimitacion-territorial-administrativa/cr_provincias_simp_wgs84.geojson",
    quiet = TRUE
  )

primates <-
  primates %>%
  st_join(cantones["canton"])

Introducción

Se van a trabajar las estadísiticas de cuatro especies de primates en Costa Rica y a su vez se mapearan los diferentes registros que se poseen.

primates %>%
  st_drop_geometry() %>%
  dplyr::select(family, species, stateProvince, locality, eventDate) %>%
  datatable(
    colnames = c("Familia", "Especies", "Provincia", "Localidad", "Fecha"),
    options = list(
      searchHighlight = TRUE,
      language = list(url = '//cdn.datatables.net/plug-ins/1.10.11/i18n/Spanish.json')
    )
  )

primates %>%
  plot_ly(
    labels = ~ etiquetas,
    values = ~ cantidades,
    type = "pie",
    textposition = "inside",
    textinfo = "label+percent"
  ) %>%
  config(locale = "es") %>%
  layout(
    title = "Especies de primates",
    xaxis = list(
      showgrid = FALSE,
      zeroline = FALSE,
      showticklabels = FALSE
    ),
    yaxis = list(
      showgrid = FALSE,
      zeroline = FALSE,
      showticklabels = FALSE
    )
  )

## Warning: 'pie' objects don't have these attributes: 'x', 'y', 'mode'
## Valid attributes include:
## 'type', 'visible', 'showlegend', 'legendgroup', 'opacity', 'name', 'uid', 'ids', 'customdata', 'meta', 'hoverlabel', 'stream', 'transforms', 'uirevision', 'labels', 'label0', 'dlabel', 'values', 'marker', 'text', 'hovertext', 'scalegroup', 'textinfo', 'hoverinfo', 'hovertemplate', 'texttemplate', 'textposition', 'textfont', 'insidetextorientation', 'insidetextfont', 'outsidetextfont', 'automargin', 'title', 'domain', 'hole', 'sort', 'direction', 'rotation', 'pull', '_deprecated', 'idssrc', 'customdatasrc', 'metasrc', 'labelssrc', 'valuessrc', 'textsrc', 'hovertextsrc', 'hoverinfosrc', 'hovertemplatesrc', 'texttemplatesrc', 'textpositionsrc', 'pullsrc', 'key', 'set', 'frame', 'transforms', '_isNestedKey', '_isSimpleKey', '_isGraticule', '_bbox'

Araña <- primates %>%
  dplyr::select(species,
                stateProvince,
                locality,
                eventDate,
                decimalLatitude,
                decimalLongitude) %>%
  filter(species == 'Ateles geoffroyi')

Carablanca <- primates %>%
  dplyr::select(species,
                stateProvince,
                locality,
                eventDate,
                decimalLatitude,
                decimalLongitude) %>%
  filter(species == 'Cebus capucinus')

Aullador <- primates %>%
  dplyr::select(species,
                stateProvince,
                locality,
                eventDate,
                decimalLatitude,
                decimalLongitude) %>%
  filter(species == 'Alouatta palliata')

Ardilla <- primates %>%
  dplyr::select(species,
                stateProvince,
                locality,
                eventDate,
                decimalLatitude,
                decimalLongitude) %>%
  filter(species == 'Saimiri oerstedii')

alt <- getData(
  "worldclim",
  var = "alt",
  res = .5,
  lon = -84,
  lat = 10
)
altitud <-
  alt %>%
  crop(provincias) %>%
  mask(provincias)

rcol <- colorNumeric(c("#9933FF", "#99FF33", "#FF3300"),
                     values(altitud),
                     na.color = "transparent")

geoffroyi <- paste0(
  "<b>",
  "Especie del primate: ",
  "</b>",
  (Araña$species),
  "<br>",
  "<b>",
  "Provincia: ",
  "</b>",
  (Araña$stateProvince),
  "<br>",
  "<b>",
  "Cantón: ",
  "</b>",
  (Araña$locality),
  "<br>",
  "<b>",
  "Fecha: ",
  "</b>",
  (Araña$eventDate),
  "<br>",
  "<b>",
  "Cordenadas: " ,
  "</b>",
  (Araña$decimalLongitude),
  "<br>",
  (Araña$decimalLatitude)
)
oerstedii <- paste0(
  "<b>",
  "Especie del primate: ",
  "</b>",
  (Ardilla$species),
  "<br>",
  "<b>",
  "Provincia: ",
  "</b>",
  (Ardilla$stateProvince),
  "<br>",
  "<b>",
  "Cantón: ",
  "</b>",
  (Ardilla$locality),
  "<br>",
  "<b>",
  "Fecha: ",
  "</b>",
  (Ardilla$eventDate),
  "<br>",
  "<b>",
  "Cordenadas: " ,
  "</b>",
  (Ardilla$decimalLongitude),
  "<br>",
  (Ardilla$decimalLatitude)
)
palliata <- paste0(
  "<b>",
  "Especie del primate: ",
  "</b>",
  (Aullador$species),
  "<br>",
  "<b>",
  "Provincia: ",
  "</b>",
  (Aullador$stateProvince),
  "<br>",
  "<b>",
  "Cantón: ",
  "</b>",
  (Aullador$locality),
  "<br>",
  "<b>",
  "Fecha: ",
  "</b>",
  (Aullador$eventDate),
  "<br>",
  "<b>",
  "Cordenadas: " ,
  "</b>",
  (Aullador$decimalLongitude),
  "<br>",
  (Aullador$decimalLatitude)
)
capucinus <- paste0(
  "<b>",
  "Especie del primate: ",
  "</b>",
  (Carablanca$species),
  "<br>",
  "<b>",
  "Provincia: ",
  "</b>",
  (Carablanca$stateProvince),
  "<br>",
  "<b>",
  "Cantón: ",
  "</b>",
  (Carablanca$locality),
  "<br>",
  "<b>",
  "Fecha: ",
  "</b>",
  (Carablanca$eventDate),
  "<br>",
  "<b>",
  "Cordenadas: " ,
  "</b>",
  (Carablanca$decimalLongitude),
  "<br>",
  (Carablanca$decimalLatitude)
)

primates %>%
  leaflet() %>%
  addProviderTiles(providers$OpenStreetMap.Mapnik, group = "OpenStreetMap") %>%
  addProviderTiles(providers$Stamen.TonerLite, group = "Stamen Toner Lite") %>%
  addProviderTiles(providers$Esri.WorldImagery, group = "Imágenes de ESRI") %>%
  addRasterImage(altitud,
                 colors = rcol,
                 opacity = 0.8,
                 group = "Altitud") %>%
  addCircleMarkers(
    data = Araña,
    stroke = F,
    radius = 4,
    fillColor = '#CC0033',
    fillOpacity = 1,
    popup = geoffroyi,
    group = "Ateles geoffroyi"
  ) %>%
  addCircleMarkers(
    data = Ardilla,
    stroke = F,
    radius = 4,
    fillColor = '#FF3399',
    fillOpacity = 1,
    popup = oerstedii,
    group = "Saimiri oerstedii"
  ) %>%
  addCircleMarkers(
    data = Aullador,
    stroke = F,
    radius = 4,
    fillColor = '#33FF66',
    fillOpacity = 1,
    popup = palliata,
    group = "Alouatta palliata"
  ) %>%
  addCircleMarkers(
    data = Carablanca,
    stroke = F,
    radius = 4,
    fillColor = '#33FFFF',
    fillOpacity = 1,
    popup = capucinus,
    group = "Cebus capucinus"
  ) %>%
  addLayersControl(
    baseGroups = c("OpenStreetMap", "Stamen Toner Lite", "Imágenes de ESRI"),
    overlayGroups = c(
      "Ateles geoffroyi",
      "Cebus capucinus",
      "Alouatta palliata",
      "Saimiri oerstedii",
      "altitud"
    )
  ) %>%
  addMiniMap(
    tiles = providers$Stamen.OpenStreetMap.Mapnik,
    position = "bottomleft",
    toggleDisplay = TRUE
  )